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Abstract 

For many physical quantities, theory supplies weak- and strong-coupling expansions 
of the types ^a n a n and a p ^2 b n (a~ 2 l q ) n , respectively. Either or both of these may 
have a zero radius of convergence. We present a simple interpolation algorithm which 
rapidly converges for an increasing number of known expansion coefEcients. The 
accuracy is illustrated by calculating the ground state energies of the anharmonic 
oscillator using only the leading large-order coefficient 60 (apart from the trivial 
expansion coefficent ao = 1/2). The errors are less than 0.5% for all g. The algorithm 
is applied to find energy and mass of the Fròhlich-Feynman polaron. Our mass is 
quite different from Feynman's variational approach. 
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1) Recently, the Feynman-Kleinert variational approximation to path integrals [1] has been 
extended to a systematic variational perturbation expansion [2]. This expansion converges 
uniformly and fast (for the anharmonic oscillator like e - constxNl/3 ' m the order N of the 
approximation [3]). Due to the uniformity of the convergence, it has given rise to an efficient 
method for extracting strong-coupling expansions from a weak-coupling expansions [4-6]. 

For many physical systems, there exists an independent knowledge of expansion coeffi- 
cients for weak and strong couplings. Important examples are most lattice models of sta- 
tistical mechanics. The purpose of this note is to propose a simple algorithm by which the 
variational perburbation expansion can be used to find a systematic convergent interpolation 
between weak- and the strong-coupling expansions. 

The algorithm is completely general and holds for any physical system whose quantities 
possess expansions in some coupling constant a of type ^a n a n , for weak and of the type 
& p J2 i> n (a~ 2 l q ) n for strong couplings, where either or both of these expansions may have a 
zero radius of convergence. 

A typical example is the ground state energy of the anharmonic oscillator with p = 
1/3, q = 3. We shall use this example to illustrated the power of the algorithm by calculating 
this the leading large-order coefficient 6 , apart from the trivial coefficent a = 1/2. The 
errors are less than 0.5% for all g. 

To make a prediction, we apply the algorithm to the Fròhlich-Feynman polaron [8], [9], 
where we know for the ground state energy the lowest three perturbation coefficients for 
weak couplings as well as the first two strong-coupling coefficients; for the polaron mass, 
only two strong-coupling coefficients have been calculated. 

Our interpolation results can be compared with Feynman's famous variational solution 
with interesting discrepàncies for the mass, calling for a calculation of higher perturbation 
or strong-coupling coefficients. 



2) Following the method explained in [3], we rewrite the weak coupling expansion of or- 
der N, 

2 



N 

En = ^2 a n° r ' ''■> 

n=0 



as 



N - a 



ra=0 



where u is an auxiliary parameter whose value is eventually set equal to 1 and p } q are two 
parameters to be determined by general properties of the strong-coupling expansion. We 
now replace u by the identical expression 



Vn 2 + u 2 -n 2 (3) 



and reexpand E N in powers of A treating u 2 — 2 as a quantity of order a. The reexpanded 
series is truncating after the order n > N. 
The resulting expansion has the form 

N / a \" 
W N (a,ü)=ü^ a Mn)( — ) (4) 

n=l 

where 

^^S( (P T ,/2 )«-»'( 1 -^· 

Forming the first and second derivatives of Wjv(«, with respect to we find the positions 
of the extrema and the turning points. The smallest among these is denoted by Ojy- The 
resulting Wn(q) = Wjv(«,^jv) constitutes the desired approximation to the energy. 

It is easy to take this approximation to the strong-coupling limit a — ï oo. For dimensional 
reasons, Ojv increases with a like Ojv ~ a^cjv, so that 

^(«^ivl^^V^J» (6) 

where 



w 



(o) 

N 



N 



J2 anfn(oo) [ — ) . (7) 



n=0 V-N, 



The full strong-coupling expression is obtained by writing Wjv(a,0) = ÍI p wn((x, lü 2 /ÍÏ 2 ), 
with à = a/0 g , and expanding w n in powers of üj 2 /Í1 2 . which behaves for a — y oo like 
(l/c 2 )(a/u q )- 2 / q . The result is 

2/g / n \ -4/ 3 



W N (a) = a p/q 
with 
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bjc) = ^wff(à,0)à^-^% =1/cq (9) 



and the superscript (n) denotes the nth derivative with respect to ú 2 . 

The parameters p and ç in the expansion (2) are now determined to render the correct 
leading and the successive powers of a in the strong-coupling expansion (8). 

The leading coefficient cjv in the optimal frequency íl n is found by searching for the 
extrema of the leading coefficient b (c) as a function of c and choosing the smallest of them. 

Explicitly 



j-n^p-lq-2n (íü) 



e 



Next we have to correct for the fact that for large but finite a has corrections to the 
behavior a l l q c. The coefficient c will depend on a like 

c(a) = C + Cl (-j +c 2 (-j +..., (11) 

requiring a reexpansion of c-dependent coefficients b c n in (8). The expansion coefficients 7„ 
are determined by extremizing b 2n (c). The final result can again be written in the form (8) 
with b c n replaced by b n which are determined by the equations shown in Table 1. The two 
leading coefficients receive no correction and are omitted. 

It is now obvious that the knowledge of any strong-coupling coefficients 6 , b\ can be ex- 
ploited to determine approximately further coefficients ajy+i, ajv+2, • • • and thus carry Wjv(a) 
to higher orders. We merely have to solve Eq. (10) for as many b n as are available. 
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3) The weak-coupling expansion of the anharmonic oscillator looks like (1) with a = g/4 
(for a potential is gx 4 /4). The lowest coefficient has is given by the ground state energy and 
has the value a = 1/2. 

The strong-coupling behavior is known from general scaling arguments to start out 
like g 1 / 3 followed by powers of g -1 / 3 , </ _1 , g~ 5 ^ 3 - Inspection of (8) shows that this corre- 
sponds to p = 1 and q = 3. The leading coefficient is known extremely accurately [7,5], 
b = 0.667 986 259 155 777108 270 962 016 919 860 .... This is now used to determine an 
approximate a\ (forgetting that we know the exact value a\ x = 3/4). The energy (4) reads 
for N = 1: 

+ 4 + < 12 > 

Equation (10) yields, for n = 0: 

bo = - + -+. (13) 

Differentiating b with respect to c we find c = c\ = 2a| / ' 3 . Inserting this into (13) fixes 
ai = (4/36o) 3 = 0.773 970 . . ., quite close to the exact value. With our approximate a\ we 
calculate Wi(a } íl) at its minimum, where 



-^u cosh [|acosh(#/#( n ))] g > g H ^ 

for (14) 

-^u cos [i arccos(^/^( n ))] 9 < 9 {n \ 



with g(°) = 2cü 3 /3v / 3· The result is shown in Fig. (1). Since the difference with respect to 
the exact solution would be to small to be visible on a direct plot of the energy, we display 
the ratio with respect to the exact energy Wi(a)/E ex . The accuracy is everyhwere better 
than 99.5 %. For comparison, we also display the much worse (although also quite good) 
variational perturbation result using the exact a ex = 3/4. 

4) Let us now turn to the polaron model. The Hamiltonian operator reads 

H = ^ + J2 ^o4«k + E {V^e lkx + h.c.) (15) 



where m& is the effective mass of the electron in the conduction band, p is the electron 
momentum, c^o is the frequency of optical phonons which are created and annihilated by a k 
and ak, and 

V4 = -'m(— ) te) (16) 

specifies the electron-phonon interaction in the volume V. The Fròhlich coupling constant 



Q = ^L ^ í±_ _ M . 17 n 

hc V 2híü \£oo £oJ 

involves the fundamental constants e } c } h and the electrònic and static dielectric constants 
£00 and e , respectively. This form of Vk assumes the size of the polaron to be large with 
respect to the lattice spacing. It further ignores spin and relativistic effects and the dispersion 
of the electron band. 

In natural units with h = c = rríf, = üjo = 1, the partition function of the polaron in 
thermal equilibrium at a fixed temperature T is described by the path integral 

1 rP , .0 a rP rP , , e -l T - T 'l 



/i [P a f> J f> J e 1 1 

Vx{T)eXP "2 L dT * 2 + ¥ïlJo Jo dTdT ' w(r)-^) 
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where fí = 1/T is the inverse temperature (at Boltzmann constant &b = 1). The weak- 
coupling expansion of the energy of the polaron is known up to the order a 3 [10]: 

E w = -a- 0.0159196220a 2 - 0.000806070048a 3 - 0(a 4 ). (19) 

For strong couplings, the energy is [11] 

E s = -0.108513a 2 - 2.836 - 0(a" 2 ). (20) 

The polaron mass has the corresponding expansions [12], [11]: 

a 

m w = 1 + - + 0.02362763a 2 + 0(a 4 ) (21) 
6 

m s = 0.0227019a 4 + 0(a 2 ). (22) 

Feynman was the first to find a uniform all-coupling constant expressions from a variational 
approximation to the path integral (18): 
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and 



the latter being evaluated at the parameters u(a), u;(a;) obtained in minimizing E F . For weak 
coupling Feynman's expressions are exact only to the order a. They have the expansions 
[10]: 

E F ' W = -a - 0.012345a 2 - 6.43434 x 10" 4 a 3 - 4.643 x 10" 5 a 4 - 3.93 x 10" 6 a 5 - . . . (25) 
m F ' w = 1 + ^ + 2.469136 x 10" 2 a 2 + 3.566719 x 10" 3 a 3 + 5.073952 x 10" 4 » 4 + . . . (26) 

For strong couplings, the expansions are 

E F ' S w -0.106103a 2 - 2.8294 - 4.86399/a 2 - 34.1952/a 4 + . . . (27) 
m F ' s w 0.020141a 4 - 1.012775a 2 + 11.85579 + . . . (28) 

With the help of the interpolation algorithm based on the variational perturbation expan- 
sion we shall find new expressions for E and m which share with Feynman's the validity for 
all a, but are more reliable at small and large a by possessing the presently most precise 
weak- and strong-coupling expansions (19), (20) and (21), (22). 

5) We now apply our interpolation algorithm the expansions (19) and (20) for the energy. 
To make the series start out with a° as required by the general ansatz (2), we remove an 
overall factor — a from E and deal with —E/a. 

Then we see from (20) that the correct leading power in the strong-coupling expansion 
requires taking p = l,ç = 1. The knowledge of b and bi allows us to extend the known 
weak coupling expansion (19) by two further expansion terms. Their coefficients a 3 ,a 4 are 
solutions of the equations 



35 15 a 2 2a 3 cia 

b = a c + (*! + - + — r + (29) 

128 8 c c 2 c 3 y J 

35 (Iq 5 ü>2 ^3 / \ 

61 = 327 " 4^ + "^- (30) 



The constant c governing the growth of Ojv for a — y oo is obtained by extremizing b in c, 
which yields the equation 

35 15 a 2 4a 3 4a 4 

Ï28 a °-Y^- — " — = °- (31) 
The simultaneous solution of (29) — (31) renders 

c 4 = 0.09819868, 

a 3 = 6.43047343 x 10" 4 , (32) 
a 4 = -8.4505836 x 10 -5 . 

The reexpanded energy (4) reads explicitly (for E including the earlier-removed factor —a) 

í 35 o 35 35 7 5 \ 2 

Wa \ol, \ l = anà it — H tttt -— H — — aia 

v ' ; V 128 320 640 3 320 5 1280 7 / 

3 / 15 5 ~3\ 4Í 2 1 \ 5 1 

+ a 2 a 3 1 H + a 3 a 4 H - a 4 a 5 — - (33) 

Extremizimg this we find 4 as a function of a [it turns out to be quite well approximated 
by the simple function 4 sa c 4 a + 1/(1 + 0.07a)]. This is to be compared with the optimal 
frequency obtained from minimizing the lower approximation W^o:, fi): 



nl = l + ^x* + J(l + ^x>) -1, (34) 
3a V V 3a 



which behaves likes c 2 a + 1 + . . . with c 2 = y8a 2 /3a ~ 0.120154. The resulting energy is 
shown in Fig. 2, where it is compared with the Feynman variational energy. For complete- 
ness, we have also plotted the weak-coupling expansion, the strond-coupling expansion, the 
lower approximation W 2 {o) , and two Padé approximants which were given in the last of 
Refs. [10] as upper and lower bounds to the energy. 

6) Consider now the polaron mass, where the strong-coupling behavior (22) fixes p = 4, q = 1. 
The coefficient b allows us to determine an approximate coefficient a 3 and to calculate the 
variational perturbation expansion W 3 (a). From (10) we find the equation 

b = —ciic 3 /8 + a 3 c, (35) 



whose minimum lies at c 3 = i/8a 2 /3ao [this value follows, of course, also directly from (37], 



where b = ^32^/27^. Using b from (22), we obtain a 3 = [27a 1 bl/32] 1 ^ 3 w 0.0416929. 
From (4) we have 

/ O 3 30 3 \ 

W 3 (a, üj) = a + dia í -— + — + — 1 + a 2 a 2 + a 3 a 3 Q. (36) 
This is extremal at 



«i = 1 + ÍÍV + J(l + ^) - 1. (37) 

3ai V V 3ai / 

The approximation W^a) = W 3 (a } íl 3 ) for the polaron mass is shown in Fig. 3, where it 
is compared with the weak and strong-coupling expansions and with Feynman's variational 
result. To see better the differences between the strongly rising curves, we have divided out 
the asymptotic behavior M as = 1 -\-b a 4 before plotting the data. As for the energy, we have 
again displayed two Padé approximants given by the last of Refs. [10] as upper and lower 
bounds to the energy. Note that our interpolation differs considerably from Feynman's and 
higher order expansion coefficients in the weak or the strong coupling expansions will be 
necessary to find out which is the true behavior of the model. 
Our curve has, incidentally, the strong-coupling expansion 

m s = 0.0227019a 4 + 0.125722a 2 + 1.15304 + 0(a" 2 ), (38) 

the a 2 -term being in sharp contrast with Feynman's expression (28). On the weak-coupling 
side, a comparison of our expansion with Feynmans's in Eq. (26) shows that our coefficent 
a 3 sa 0.0416929 is about 10 times larger than his. 

Both differences are the reason for our curve forming a positive arch in Fig. 2, whereas 
Feynman's has a valley. It will be interesting to find out how the polaron mass really 
behaves. This would be possible by calculating a few more terms in either the weak- or the 
strong-coupling expansion. 

Note that our interpolation algorithm is much more powerful than Padé's. First, we 
can account for an arbitrary fractional leading power behavior a p as a — ï 00. Second, the 
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successive lower powers in the strong-coupling expansion can be spaced by an arbitrary 2/q. 
Third, our functions have in general a cut in the complex a-plane approximating the cuts 
in the function to be interpolated [13]. Padé approximants, in contrast, have always an 
integer power behavior in the strong-coupling limit, a unit spacing in the strong-coupling 
expansion, and poles to approximate cuts. 
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FIGURES 




FIG. 1. Plot of the ratio of the intepolation energy with respect to the exact energy as a function 
of the coupling constant. The accuracy is everyhwere better than 99.5 %. For comparison, we also 
plot the variational perturbation result using the exact a\ = 3/4. 
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FIG. 2. The polaron energy obtained from our variational interpolation between the weak 
coupling expansion (dashed) and the strong-coupling expansion (short-dashed) to be compared 
with Feynman's variational approximation (fat dots). The dotted curves are Padé approximants 
which were given in the last of Refs. [10] as upper and lower bounds to the energy. 
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FIG. 3. The polaron mass curve interpolating optimally between the weak-(dashed) and 
strong-coupling expansions (short-dashed). To see better the differences between the strongly 
rising functions, we have divided out the asymptotic behavior M as = 1 + boa 4 before plotting 
the curves. The fat dots show Feynman's variational approximation. The dotted curves are Padé 
approximants which were given in the last of Refs. [10] as upper and lower bounds to the mass. 



13 



TABLES 

TABLE I. Equations determining the coefficients b n in the strong-coupling expansion from the 
functions b n (c) and their derivatives. For brevity, we have suppressed the argument c in the entries. 



n 


K 


-c n -\ 


2 
3 

4 


b 2 + ci&i + ±cX 

b 3 + c 2 b[ + ci6' 2 + Cl c 2 b% + ic?6ï + ±c?&( 3) 
& 4 + c 3 6i + c 2 b' 2 + ci 6^ + (|c| + c lC3 )b% 
+ Cl c 2 b'{ + i C 2ò' 2 ' + lclc 2 b^ + icfòf ) + ±cibP 


b' 2 + Cí b'{+\clb^)/b^ 
{b'z + c^'i + c^ + c^bf 
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